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The monopole order parameter of QCD is computed in terms of gauge invariant field strength correlators. Both 
quantities are partially known from numerical simulations on the lattice. A new insight results on the structure 
of the confining vacuum. 



1. Introduction 

The mechanism of confinement by dual su- 
percoductivity of the QCD vacuum PP|21I3| is 
confirmed by numerical simulations on a lattice. 
Chromoelectric flux tubes connecting q — q pairs 
produced by dual Meissner effect are indeed ob- 
served in lattice configurations with the expected 
form of the electric and magnetic fields @]. An 
extensive analysis has been performed by explor- 
ing the vacuum by means of an order parame- 
ter (fi) |5l6l7l8l5] which is the vacuum expecta- 
tion value of a magnetically charged operator /i. 
In the confined phase (/i) ^ , which implies 
Higgs breaking of the magnetic U(l) symmetry, 
in the deconfined phase (fi) = and the magnetic 
charge is superselected. In SU(N) gauge theory 
there exist — 1 magnetic charges and A^ — 1 in- 
dependent operators n"" , {a — 1, .., N — 1), which 
create monopoles of the species a PHI- They can 
be written 

q is an integer b± is the field of a Dirac monopole 
with — and V A 6x = -3- + Dirac — string. 



(2) 



with U{x, y) an arbitrary gauge transformation. 
We will take for U a parallel transport to x from 
a reference point y along a path C . 
(/i") is gauge invariant. In Eq(2) 



*3 = diag{ 



< a > 

A-a N-a 



N 



N 



< -(N-a) > 



In the gauge in which $° = <I>J^ (Abelian Projec- 
tion) 



= / d^yb^{x-y)Eliy.t) 



(4) 



where E'^ is the transverse chromoelectric field 
along the color direction T° 

a . a + I 
^ diag{0,...0,l , -1,0, ..0) (5) 

E'^ ± is the conjugate momentum to A'^j_. Hence 

^^'^ix■,t)\Aliy,t)) = \Aliy,t) + ^b^{x ^ y)) (6) 

/i° {x, t) creates a Dirac monopole at {x, t) in the 
residual gauge symmetry after abelian projection. 
It proves convenient to use instead of (/z") the 
susceptibility 



(7) 



Here (3 is the usual variable of the lattice formu- 
lation f3=^. Since at /3 = = 1, 



(8) 



has been measured by lattice simulation for 
various gauge theoriesxompact U{1) ^,SU{2) [S] 
, SU{3) and Nf ^ 2 QCD |TT] . In all these 
systems — > finite in the confined phase in the 
thermodynamical limit V = 00. By use of 

Eq(8) this means (//°) ^ . 
In the deconfined phase T > Tc 



« -\c\Ls + c' or {fi'') = 



(9) 
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or, again by Eq(8) {n'^) = . 

In the critical region T k, T^. the scahng law 
holds 



o 1 



(10) 



Here r = 1 — -^ , v is the critical index of the cor- 
relation length of the order parameter, is in- 
dependent of the choice of the abelian projection 
[T^ |14j . Expanding the exponential which 
defines one has 

- y„)(($'^J).„,j;,($«.^).„,,j;„) (11) 

The notation is {^".E),^^ = Tr[^°'{x).E^{x). 

The vev's in Eq(ll) are gauge invariant field 
strength correlators. 

We shall identify these correlators with those 
of the Stochastic Vacuum approach to QCD [501 
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2. Cluster Expansion of (/i") 

In QCD all observables can be expressed in 
terms of gauge invariant field strength correlators. 
The basic idea of the stochastic vacuum approach 
is to perform a systematic cluster expansion of 
the correlators and truncate it typically retaining 
only the clusters up to order two. Since the one 
point cluster (($'^.£')) — 0, only two point clus- 
ters will survive the truncation, and hence only 
even terms in the expansion Eq(ll) , which will 
be approximated as products of two point corre- 

lators ^i^.^im-m) = ((1>^^)n,yl($^^).„y.). 

There is a combinatorial factor (2n — 1)!! for the 
term of order 2n . Eq(ll) becomes 



These correlators are measured on the lattice 
QTldlCSI I and are used as an input in stochas- 
tic QCD. For those correlators {E^'E'') = 5'^''$ , 
so that p° is independent on a . This is also the 
case in lattice determinations of p° 0. 

The the cluster expansion is generically ex- 
pected to work at large distances, and in the 
study of confinement we are looking for infrared 
properties . Anyhow a direct check of it can be 
obtained by looking at the dependence of on 
q. The truncated is proportional to : higher 
correlators would introduce terms proportional to 
higher powers of g. Old data [H] 53 seem to agree 
with but a systematic study of this dependence 
will be done. 

3. The Field Correlators 

A general paranietrization of field strength cor- 
relators dictated by invariance arguments [20! [21] 
is 

*f „ (zi - Z2) = — (TrF" „ (zi) 

V{zi,Z2)Ff,,,,,iz2)V\zi,Z2)) (14) 

*f . „ (21 - Z2) = 5"'' 

1 d 

+ 7^-^ [Di{zi - Z2){Z^,J^^^., - Z^^J^^f,^] + 



1 d 



2dz 



[D{Z,-Z2){ 



At r ^ the electric field correlators do not co- 
incide with the magnetic ones, and there are four 
form factors, De, Die, Djj, Dm ■ 

For correlators of electric fields Ei^ , Ei^ Eq(15) 
gives 



1 



(p-) = e-^^/'i^w/d-'y2*?,„,fo-i-y2)b-(yi)fc-(y2)^^2) *M2 = S-%i;2{De + ^D^) + 



d 
dii. 



(16) 



or, since (3 



2N 



P = 



q^ d 



[/3 / d^yi / dPy2 



(13) 



In the convolution with b± the derivative terms 
give 0. For the same reason 

^1-112 ^ ^2l«2 fc^ 

Going to the Fourier transform we get for p° 
Eq(13) 



We shall identify with the two point correla- 
tors defined with a straight line parallel transport. 



lQdl3^' J (27r)3 
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DEik')^{k^d,,,,-kM] (17) 
Here DE{k^) is the Fourier transform of {De + 
Since 

{^5,, - hih2)blik)b'l{~k) = |i?(fc)|2 (18) 

we can use the exphcit form of H{k) 

H{k) = kAb^{k) (19) 

with n the direction of the Dirac string ( we shaU 
cah it 2;), and get 



\H{k)\' 



1 1 



(20) 



For we then have 

where /(t^) s ^De(I!^-). 

Identifying our correlators with those of the 
stochastic model simply explains that p° is in- 
dependent both on a and on the abelian projec- 
tion. At large /3(deconfined phase) /(fc^) can be 
approximated by first order perturbation theory 



fik') - 



2Nk 



(22) 



The only dependence on /3 is the explicit factor 
in Eq(21) so that 

^ " 167V 7 (2^)3 fc^P (^'^) 

The integral is easily computed with UV cut-off 
i (a the lattice spacing) and IR cut-off [Lg 
the spacial size of the lattice) . The result is 



P 



im (27r) 



^[-v/(2)L,. + 2ln{Ls) + const]{2A) 



By comparison with Eq(8) this means that (/i") = 
in the thermodynamical limit Ls —>■ 00. The 
correlators have been measured on the lattice 
both at T = and at T 7^ 0. In the range 
of distances .Ifm < x < Ifm they are well 
parametrized by a form jl7| 



De = Ae ^t. + —e' 



(25) 



with Xa ~ 2Af, and At, « .3/m A and b are in- 
dependent on P within statistical errors in the 
range from T = up to T « .95Tc . Approaching 
further Tc A rapidly decreases to zero.|19|. 

In fact the parametrization Eq(25) cannot be 
valid at shorter distances, where the operator 
product expansion and the non-existence of con- 
densates of dimension less than 4 require that 

The prescription on the singularity is the same 
as in perturbation theory. At larger distances a 
stronger infrared cut-off at some distance A, must 
exist, since colored particles cannot propagate at 
infinite distance. This feature needs a further nu- 
merical investigation of correlators at large dis- 
tance on the lattice. 

Up to T « .95Tc the only dependence on (3 in 
Eq(21) is again the explicit factor so that the re- 
sult coincides with Eq(24) except that the lattice 
size Ls is replaced by the infrared cut-off A 



1 



P 



'y/{2)A + 2ln{A) + const] (27) 



16 (27r)2 ' 

The integral of the exponential term of Eq(25) is 
included in the constant. This expression gives 
a finite value of for T < Tc independent of 
the volume and hence (/i") ^ .This means dual 
superconductivity for any finite value of the UV 
cut-off a. However in the continuum limit a — > 
p° diverges so that needs a renormalization. 
This is similar to what happens for the Polyakov 
line P2 ■ Existing Lattice data support this state- 
ment [See Fig(2) of ref.|Hl] i^ut we plan to do a 
more systematic investigation of this issue. By 
approaching the critical temperature both the IR 
cut-off A and the coefficient A in Eq(25) strongly 
depend on /3 :A diverges and A tends to zero. A 
more detailed calculation, which will be reported 
elsewhere, gives 



P 



d 



{P[j^{-Vl2)- + 2\n{-) 
a a 



16N df3^ '(27r)2 
+const) - 



2NXtAE{l + ^] (28) 



Numerical determinations of the dependence of 
Ae on of the temperature around Tc exist in the 
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literature^^. Not much is known about the be- 
havior of A . Further study is needed to under- 
stand how the scahng law of p° Eq(lO) described 
in Sect 1 conies out of Eq(28). 

4. Conclusions 

We close with a few remarks. 

Checking the dependence p° c>c is a test of 
the Stochastic Vacuum model. 

The independence of p° on a and on the abelian 
projection are also an important test of it. 

The existence of confinement depends on a 
strong infrared cut-off of the field correlators. 
This can be directly checked on lattice.. 

p° diverges in the continuum limit, but provides 
a good description of confinement at any fixed 
value of the UV cut-off. 

An interesting interplay emerges of confine- 
ment with infrared properties of gauge invariant 
field strength correlators. 

Useful discussions with H.G. Dosch, Y.A. Si- 
monov,M.D'Elia, E. Meggiolaro, G.Paffuti are 
aknowledged. 
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